Additivity of effective quadrupole moments and angular momentum alignments in the 

y4 ~ 130 nuclei 
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The additivity principle of the extreme sheU model stipulates that an average value of a one-body 
operator be equal to the sum of the core contribution and effective contributions of valence (particle 
or hole) nucleons. For quadrupole moment and angular momentum operators, we test this principle 
for highly and superdeformed rotational bands in the A ~ 130 nuclei. Calculations are done in the 
self-consistent cranked non-relativistic Hartree-Fock and relativistic Hartree mean-field approaches. 
Results indicate that the additivity principle is a valid concept that justifies the use of an extreme 
single-particle model in an unpaired regime typical of high angular momenta. 

PACS numbers: 21.60.Jz, 21.60.Cs, 21.10.Gv, 21.10.Ky, 23.20.Js, 27.60.+j 



I. INTRODUCTION 

The behavior of the nucleus at high angular momenta 
is strongly affected by the single-particle (s.p.) structure, 
i.e., shell effects. Properties of the s.p. orbits around 
the Fermi level determine the deformability of the nu- 
cleus, the amount of angular momentum available in the 
lowest-energy configurations, the moment of inertia, and 
the Coriolis coupling. Consequently, nuclconic shells can 
be seen and probed through the measured properties of 
rapidly rotating nuclei. 

The independent particle model is a first approxima- 
tion to the nuclear motion. Here, the nucleons are as- 
sumed to move independently of each other in an average 
field generated by other nucleons. Each nucleon occu- 
pies a s.p. energy level, and levels with similar energies 
are bunched together into shells. The wave function of 
a given many-body configuration uniquely characterized 
by s.p. occupations is an antisymmetrized product of one- 
particle orbitals (the Slater determinant). In the next 
step, the residual interaction between particles needs to 
be considered. This is the essence of the configuration 
interaction method or the interacting shell model. For 
heavier nuclei, where the number of s.p. orbits becomes 
large, a customary approximation is to divide the config- 
uration space into the (inert) core states and the (active) 
valence orbits and to perform configuration mixing in the 
valence subspace. 

The basic idea behind the additivity principle for one- 
body operators is rooted in the independent particle 
model. The principle states that the average value of a 



one-body operator O in a given many-body configuration 
k, 0{k), relative to the average value in the core configu- 
ration 0'^°'^'^ , is equal to the sum of effective contributions 
of particle and hole states by which the A:-th configura- 
tion differs from that of the core. Such a property is 
trivially valid in the independent particle model. How- 
ever, the presence of residual interactions and resulting 
configuration mixing could, in principle, spoil the simple 
picture. In particular, in the interacting shell model, the 
polarization effects due to additions of particles or holes 
are significant and they give rise to strong modifications 
of the mean field. So the essence of the additivity prin- 
ciple lies in the fact that these polarizations are, to a 
large extent, independent of one another and thus can 
by treated additively. 

The additivity principle for strongly deformed nuclear 
systems was emerging gradually in the 1990s. First, it 
was found in Ref. [I] that effective (relative) angular mo- 
mentum alignments are additive to a good precision in 
the superdeformed (SD) bands around ^''^Gd. However, 
the analysis was only restricted to a few bands. Later, 
the statistical analysis of Ref. [2] in the A ~ 150 and 190 
mass regions clearly demonstrated that the so-called phe- 
nomenon of band twinning (or identical bands) is more 
likely to occur in SD than in normal-deformed bands. It 
was shown that a necessary condition for the occurrence 
of identical bands is the presence of the same number of 
high-iV intruder orbitals (see also Ref. [3]). In addition, it 
was concluded that the configuration-mixing interactions 
such as pairing and the coupling to the low-lying collec- 
tive vibrational degrees of freedom act destructively on 
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identical bands by smearing out the individuality of each 
s.p. orbital. Such individuality is an important ingredi- 
ent for the additivity principle: it is expected that this 
principle works only in the systems with weak residual 
interaction, in particular, pairing [2, 4]. 

The principle of additivity at superdeformation was 
explicitly and thoroughly formulated for the case of the 
Q20 quadrupole moments in the non-relativistic study of 
quadrupole moments of SD bands in the A ~ 150 mass 
region in Ref. [5] within the cranked Hartree-Fock (CHF) 
approach based on Skyrme forces. It was shown that 
the charge quadrupole moments calculated with respect 
to the doubly magic SD core of ^^^Dy can be expressed 
very precisely in terms of effective contributions from the 
individual hole and particle orbitals, independently of the 
intrinsic configuration and of the combination of proton 
and neutron numbers. 

Following this work, it was shown that the princi- 
ple of additivity of quadrupole moments works also in 
the framework of the microscopic -l-macroscopic method 
(in particular, the configuration-dependent cranked Nils- 
son-1-Strutinsky approach) [6, 7]. However, contrary to 
sc;lf-c,onsistent approaches, the effective s.p. quadrupole 
moments of the microscopic-|-macroscopic method are 
not uniquely defined due to the lack of self-consistency 
between the microscopic and macroscopic contributions. 

The study of additivity of quadrupole moments and ef- 
fective alignments was also performed in the framework 
of the cranked relativistic mean field (CRMF) approach, 
but it was restricted to a few configurations in the vicin- 
ity of the doubly magic SD core of ^^^Dy [8]. It was 
suggested in this work that the additivity principle when 
applied to the angular momentum operator (i.e., effective 
afignments) does not work as well as for the riuadrupole 
moment. In addition, the effective alignments of high-A?^ 
intruder orbitals seem to be less additive than the effec- 
tive alignments of non-intruder orbitals. The latter can 
be attributed to a pronounced polarization of the nucleus 
by high- A/' intruder orbitals at high spin. 

For quadrupole moments, the additivity principle was 
experimentally confirmed in the A~140-150 mass region 
of superdeformation. It was shown that the quadrupole 
moments of the SD bands in ^^"^Sm [9] and ^'^^Gd [10] 
could be well explained in terms of the ^^^Dy SD core 
and effective s.p. quadrupole moments of valence (parti- 
cle and hole) orbits. All of these studies, together with 
the previous results for moments of inertia [11, 12] and 
effective alignments [8, 13], strongly suggest that the SD 
bands in the ^^^140-150 mass region are excellent ex- 
amples of an almost undisturbed s.p. motion. This is 
especially true at rotational frequencies above fi(x-'=G.5 
MeV [8, 12] where pairing is expected to be of minor 
importance. (For other excellent examples of an almost 
undisturbed s.p. motion at high spins, see Refs. [14 16].) 

In the mass A~135 (Z=58-62) light rare-earth re- 
gion, large Z=58 and N=72 shell gaps (see Fig. 1 and 
Refs. [17, 18]) lead to the existence of rotational struc- 
tures with characteristics typical of highly deformed and 



SD bands. These bands were observed up to high and 
very high spins (see Refs. [19, 20] and references quoted 
therein). For example, the yrast SD band in ^''^Cc ex- 
tends to ~68/i, which represents one of the highest spin 
states ever observed in atomic nuclei [19]. At such high 
spins, pairing is expected to play a minor role [11, 18, 21], 
which is a necessary condition for the additivity principle 
to hold. In this mass region, experimental studies of the 
additivity principle were performed in Refs. [20, 21]. Dif- 
ferential lifetime measurements, free from common sys- 
tematic errors, were performed for over 15 different nuclei 
(various isotopes of Ce, Pr, Nd, Pm, and Sm) at high spin 
within a single experiment [20, 21]. 

There are several notable differences between the 
A~135 and A~140-150 regions of superdeformation. In 
particular, the rotational bands in the ^~135 region are 
calculated to correspond to the local energy minima that 
are characterized by much larger 7-softness than those 
in the ^~140-150 mass region [17, 18]. Thus, one of the 
main goals of the present manuscript is to find the impact 
of the 7-softness on the additivity principle. The second 
goal is a detailed study of the additivity principle not only 
for quadrupole moments but also for angular momentum 
alignments. The present work is the first study where the 
additivity of relative alignments has been tested within 
the CHF and CRMF frameworks in a systematic way 
along with the additivity of quadrupole moments. Some 
results of this study have been reported in Refs. [20, 21]. 

This paper is organized as follows. The principle of 
additivity, definitions of physical observables, the way of 
finding effective s.p. quantities, and details of theoreti- 
cal calculations are discussed in Sec. II. Analysis of the 
additivity principle for quadrupole moments and relative 
alignments, and the discussion of associated theoretical 
uncertainties are presented in Sec. HI. Finally, Sec. IV 
contains the main conclusions of our work. 

II. THEORETICAL FRAMEWORK 

A. Definition of observables 

Since pairing is neglected in this work, the charge 
quadrupole moments Q2Q and Q22 are defined microscop- 
ically as sums of expectation values of the s.p. quadrupole 
moment operators 920 and 922 of the occupied proton 
states, i.e., 

Q20 = ^(M|92o|Ai), (1) 
ft 

Q22 = ^(m|9221m), (2) 
ft 

where 920 and 922 are defined in three-dimensional Carte- 
sian coordinates as [22] (conserved signature symmetry 
is assumed) 

520 = 2z^-x^-y\ (3) 
922 = V3(a;2-2/2). (4) 
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neutrons protons 




proton quadrupole moment (eh) 

FIG. 1: (color online) Single-particle energies for neutrons (left) and protons (right) in ^'^^Ba as a function of the proton 
quadrupole moment calculated in the HF-|-SLy4 model. Solid and dashed lines mark positive and negative parity states, 
respectively. The orbitals are labeled by the asymptotic (Nilsson) quantum numbers [f\fnzA]Q of the dominant component of 
the s.p. wave function. The neutron intruder orbitals originating from the J\f=6 shell are additionally labeled by the main 
oscillator quantum number and a subscript denoting the position of the orbital within the Af shell. 



The factor of ^/3 is included in the definition of 522 
in order to have the following expressions for the total 
quadrupole moment Q2 and the associated Bohr angle 7: 

Q2 = ^JQ^, + Qh, (5) 

tan(7) = Q22/Q20- (6) 

Note that the sums in Eqs. (1-2) run only over pro- 
ton states. The neutrons, having zero electric charge, 
do not appear in the sums explicitly, but they influ- 
ence the charge quadrupole moments indirectly via the 
quadrupole polarization (deformation changes) induced 



by occupying/emptying single-neutron states. 

It should be noted that with the definitions (1-6), the 
spherical components of the quadrupole tensor are Q20 
and Q22/V^- This fact is important for the definition 
of the so-called transition quadrupole moment Qt [23, 
24]. This moment gives the measure of the transition 
strength of the A/=2 (stretched) E2 radiation in the 
limit of large deformation and angular momentum, and it 
is proportional to the component Q22/ of the spherical 
quadrupole tensor when the quantization axis coincides 
with the vector of rotational velocity u, i.e.. 



Q20 = Dio{r,e",r)Q2o + [Di2ir,o^,r) + Di_2{r,o'^,r)]^, (7) 
^ = Di.ir, o-,r)Q2o + [Di2{r,0\ r) + Di^,{r, 0^, r)] (s) 



Here, symbols Z?^^ denote the Wigner functions [25], with their arguments ip'^ , O'^ , (j>'^ being the Euler angles 
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FIG. 2: (color online) Neutron s.p. energies (routhians) in 
the self-consistent rotating potential (CHF+SLy4) as a func- 
tion of rotational frequency for the core configuration (the 
lowest highly deformed band of ^'^^Ce). Occupied (empty) 
states are denoted by thick (thin) lines. Dotted, solid, dot- 
dashed, and dashed lines indicate parity and signature quan- 
tum numbers (tt, r) = (+, +i), { — and ( — , — i), 
respectively. The orbitals are also labeled by the asymptotic 
quantum numbers [A/'n^AJSl of the dominant harmonic os- 
cillator component of the wave function. The neutron in- 
truder orbitals originating from the A/'=6 shell are addition- 
ally marked. At intermediate rotational frequencies, the low- 
est intruder level 6i becomes occupied and this leads to the 
presence of the large gap in the spectrum at N—73. 

that rotate the z axis (the standard quantization axis 
for spherical tensors) onto the direction of the angular 
velocity. 

For the cranking axis coinciding with the y-axis of the 
intrinsic system, as is the case for the code HFODD [26, 
27] used in the present study, the Euler angles are 1^=^^ 
6'=7r/2, and (/)=7r/2, which gives: 

Q20 - -2^20-^2-^, (9) 




0.0 0.2 0.4 0.6 0.8 
rotational frequency (MeV) 

FIG. 3: (color online) Similar to Fig. 2 except for proton s.p. 
states. The large proton Z=58 gap in the s.p. spectrum is 
present at all frequencies considered. 

The second of these equations gives the definition of the 
transition quadrupole moment used in this work: 

In order to provide a link to studies that employ the 
X-axis cranking, like, e.g., Rcfs. [23, 24] and our earlier 
papers [20, 21], we repeat derivations for the Euler angles 
■0=71/2, 0=7r/2, (/)=7r, which rotate the z axis onto the x 
axis: 

^22^ f^n 

^ ^ "V8^^''"2Vf' ^''^ 

hence 

=-V3^ = '^^" + V37f ^''^ 

Although definitions (11) and (14) differ by signs of the 

second terms, values of Q""^ and Qj^"^ obtained in self- 
consistent calculations must be identical because they 
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cannot depend on the direction of the cranking axis. It 
means that values of Q22 obtained in cranking calcula- 
tions along the y and x axes have opposite signs. In what 
follows, we employ definition (11) of the transition mo- 
ment and drop the superscripts that denote the direction 
of the cranking axis, e.g.. 



(15) 
(16) 



Finally, the expectation value of the total angular mo- 
mentum J (its projection on the cranking axis) is defined 
as a sum of the expectation values of the s.p. angular mo- 
mentum operators jy of the occupied states 



Qt = Q20 ~ A/ 3 '322, 




J = (Jy) =^\fJ'\jy\^J)■ 



{17) 



The value of J can be expressed in terms of the total spin 
/ via the cranking relation [28] 



j = v/T(7TT)«/ + J. 



(18) 



B. Additivity of effective s.p. observables 

For each fc-configuration defined by occupying a given 
set of s.p. orbitals and represented by a product state 
|fc), we determine the average value 0{k) = {k\0\k) of 
a s.p. operator O. We may now designate one of these 
configurations as a reference, or a core configuration, and 
determine the relative change dO{k) = 0{k) — Q'^°'^'^ of 
the physical observable in the fc-th configuration with re- 
spect to that in the core configuration. The additivity 
principle stipulates that all these differences can be ex- 
pressed as sums of individual effective contributions 
coming from s.p. states (enumerated by index a), i.e.. 



In this way, one can label the fc-th configuration with the 
set of coefficients c{k) = {ca{k), a = 1, . . . , m}, where m 
denotes the size of s.p. space considered. The values of 
can be calculated by proceeding step by step from the 
core configuration to the configurations differing by one 
particle or one hole, tlien to the configurations differing 
by two particles, two holes, or a particle and a hole, and 
so forth, until the data set is generated which is statis- 
tically large enough to provide appreciable precision for 
o^. Had the additivity principle been obeyed exactly, 
calculations limited to one-particle and one-hole configu- 
rations would have sufficed. Since our goal is not only to 
determine values of o'^ but actually prove that the addi- 
tivity principle holds up to a given accuracy, we have to 
consider a large set of configurations and determine the 
best values of together with their error bars. 

In what follows, we consider relative changes in the av- 
erage quadrupole moments SQ2oik) and 5(522 (fc), transi- 
tion quadrupole moments SQt{k), and total angular mo- 
menta SJ{k) (sec Sec. II A), which are related to the 
effective one-body expectation values via the additivity 
principle. 

The addition of particle or hole in a specific single- 
particle orbital a gives rise to a polarization of the sys- 
tem, so the effective s.p. values, 0°^, depend not only 
on the bare s.p. expectation values, o^^'^^ = (6)^, but 
also contain polarization contributions. For example, 
the effective s.p. charge quadrupole moment glo q can be 
represented as the sums of bare s.p. charge quadrupole 
moments a = {^20) a ^^'^ polarization contributions 
„poi . 



off _ bare , pol 
'ho.a ^ 'ho.a ^ 'l2Q.a 



(20) 



Therefore, for neutron orbitals, which have vanishing 
bare charge quadrupole moments, 920,an ~ ^' 
fective charge quadrupole moments are solely given by 
polarization terms: 



eff 
y20,a„ 



„Pol 
y20,Q 



(21) 



Oik) ~ Q<=°^^ = SO{k) 



cir 



(19) 



C. Determination of effective s.p. observables 



Coefficients Cq in Eq. (19) define the s.p. content of the 
configuration k with respect to the core configuration. 
Namely, 

(i) Ca{k) = if the state a is not occupied in either of 
these two configurations, or is occupied in both of 
them, 



Once the averages of physical observables 0(k) for the 
set of Nc calculated configurations (fc = 1, . . . ,Nc) are 
determined, the effective s.p. contributions (19) are 
found by means of a multivariate least-square fit (see, 
e.g., Refs. [29, 30]). This is done by minimizing the func- 
tion of defined by 



(ii) Ca{k) = 1 if a has a particle character (it is occu- 
pied in the fc-th configuration and is not occupied 
in the core configuration), 

(iii) c„(fc) = — 1 if the state a has a hole character (it is 
not occupied in the fc-th configuration and is occu- 
pied in the core configuration). 



F [of] 



50{k) 



ofc„(fc) 



(22) 



Note that the problem is only meaningful when the num- 
ber of configurations is sufficiently large, Nc > m. Fol- 
lowing the general concept of the least-square method. 
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the partial differentiation with respect to the variables 
o'^ yields 

= ^^C„(fc)Ca-(fcKj-^(50(fc)Ca(fc) = 
a' k k 

= {Bo''^-a)o., (23) 

where a„ = '^k^^i^)'^^^^) ~ and B = WB^^'W = 

\\Ylik'^oiik)ca'(k)\\ = c^c. Solving this equation by in- 
verting the non-singular matrix B gives the solution to 
the multivariate regression problem: 

o^ff = B-^a = {c'^c)-^c'^SO. (24) 

The fact that B is positive-definite guarantees that the 
solution o'^^ corresponds to a minimum of F [o*^] . 

In order to estimate the variance, we assume that the 
first statistical moments of residuals, 

eo{k) = SO{k) -Y,Ca{k)df, (25) 

a 

are zero for all k = l,...,Nc. Consequently, d'^^ can be 
considered an unbiased estimate of o'^^. Furthermore, 
under the assumption that 

var(eo(fc)) = cr^ (26) 

for all = 1, . . . , Nc, and 

cov(eo(fc),eo(fc')) =0 (27) 

for all {k,k' = l,...,Nc\k ^ A:'}, one can define the 
variance-covariance matrix as a^B~^ = (7^(c"^c)~^, for 
which the unbiased estimate for cr^ is given by 

° k=l 

Finally, the unbiased estimate for the variance-covariance 
matrix for o^^ is given by -B"^^^. In what follows we 

do not differentiate between notations for variables and 
their estimates. The least-square procedure described 
in this section was used to determine the effective s.p. 
quadrupole moments {<720,a! 'i'22.ai = and 

angular momentum alignments {j^, a = 1, m}. 



D. Method of calculations 

The CHF calculations were performed using the code 
HFODD (vl.75) [26, 27] with the interaction SLy4 [31]. 
The accuracy of the harmonic oscillator (HO) expansion 

depends on the frequencies {hiUx, hujy, and fv-Oz) of the 
oscillator wave functions and the number M of the HO 



states included in the basis. The basis set includes the 
lowest M states with energies given by 

^n^,ny,n, = fu^x{nx + \) + f^yi^y + ^) + f^z{nz + 

(29) 

An axially symmetric basis {wx = cOy) with the defor- 
mation q = oJx/lUz = 1.81, oscillator frequency hujo = 
41^-1/3 MeV, and value of M = 296 was used in ah 
the CHF calculations. This basis provides sufficient nu- 
merical accuracy for the physical observables of interest 
[32]. 

The CRMF calculations were performed using the 
computer code developed in Refs. [12, 33, 34]. An 
anisotropic three-dimensional harmonic oscillator basis 
with deformation (/3o — 0.4, 7 — 0°) has been used in 
the CRMF calculations. All fermionic states below the 
energy cutoff < 11.5ftw,f and all bosonic states 

below the energy cutoff < 16.5hujQ were used 

in the diagonalization and the matrix inversion. This 
basis provides sufficient numerical accuracy. The NLl 
parametrization of the RMF Lagrangian [35] has been 
used in the CRMF calculations. As follows from our pre- 
vious studies, this parametrization provides reasonable 
s.p. energies for nuclei around the valley of /3-stability 
[8, 36]. 

E. Selection of independent-particle configurations 

In both CHF and CRMF calculations, the set of 
independent-particle configurations in nuclei around 
-'^^-'^Ce was considered. The final sets used in additiv- 
ity analysis consisted of 183 and 105 configurations in 
the CHF and CRMF variants, respectively. All ambigu- 
ous cases, due to crossings, convergence difficulties, etc., 
were removed from those sets. Since the CRMF calcula- 
tions are more time-consuming than the CHF ones, the 
CRMF set is smaller. Nonetheless, the adopted CRMF 
set is sufficiently large to provide reliable results. To put 
things in perspective, in Rcf. [5], where the CHF analysis 
of additivity principle in the SD bands of the A ~ 150 
mass region was carried out, 74 calculated SD configura- 
tions were considered. 

Every calculated product-state configuration was la- 
beled using the standard notation in terms of parity- 
signature blocks [7V+^+i, A^+,_i, iV__+i, A^_,_i], where 
NT^^r are the numbers of occupied s.p. orbitals having par- 
ity TT and signature r. In addition, the s.p. states were 
labeled by the Nilsson quantum numbers and signature 
[AfnzA]fl'^ of the active orbitals at zero frequency. The 
orbital identification is relatively straightforward when 
the s.p. levels do not cross (or cross with a small interac- 
tion matrix element), but it can become ambiguous when 
the crossings with strong mixing occur. In some cases, it 
was necessary to construct diabatic routhians by remov- 
ing weak interaction at crossing points. Even with these 
precautions, a reliable configuration assignment was not 
always possible; the exceptional cases were excluded from 
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the additivity analysis. Clearly, the likelihood of the oc- 
currence of level crossings is reduced when the s.p. level 
density is small, e.g., in the vicinity of large shell gaps. 

Large deformed energy gaps develop at high rotational 
velocity for Z=58 and N=73 (see Figs. 2 and 3). There- 
fore, the lowest SD band {i'i-\_3/2 band) in ^^^Ce is a nat- 
ural choice for the highly deformed core configuration in 
the >1~130 mass region. The additivity analysis was per- 
formed at a large rotational frequency of ?iw=0.65MeV. 
This choice was dictated by the fact that (i) at this fre- 
quency the pairing is already considerably quenched, and 
(ii) no level crossings appear in the core configuration 
around this frequency (cf. Figs. 2 and 3). Moreover, at 
this frequency, the lowest neutron ^13/2 intruder orbital 
already appears below the N^73 neutron gap (see Fig. 
2). The choice of an odd-even core, strongly motivated by 
its doubly closed character at large deformations/spins, 
does not impact the additivity scheme, which is insensi- 
tive to the selection of the reference system. 

The highly deformed core configuration in ^^^Ce 
([18,19,18,18]„0 [14,14,15,15]p) has the following orbital 
structure: 

|core) = |core)j^ (g) |core)^ = 



(^.(1*13/2)6^) 

(z.(l/in/2)[523]7/2±')2 
(Kl5i/2)[411]l/2±«)2 
(Kl.97/2)[413]5/2±0' 
(K2/,/2)[541]l/2±0' 
(1.(24/2) [411]3/2±0' 

(j.(l/^n/2)[532]5/2±0' 
(1.(159/2) [404]9/2±0' 
( ••• )|0). 



(7r(l99/2)[404]9/2± 
(^(l/*n/2)[541]3/2± 

(^(1,97/2) [420] l/2± 

(^(24/2)[422]3/2± 
(7r(l/in/2)[550]l/2±; 

(7r(2pi/2)[301]l/2±0 
( ••• )|0). 



where dots denote the deeply bound states and \0)^ and 
|0)^ are the neutron and proton vacua, respectively. The 
spherical subshells from which the deformed s.p. orbitals 
emerge (cf. Fig. 1) are indicated in the front of the Nilsson 
labels. 

The Nilsson orbital content of an excited configuration 

is given in terms of particle and hole excitations with 
respect to the core configuration through the action of 
particle/hole operators with quantum labels correspond- 
ing to the occupied or emptied Nilsson orbitals. The 
character of the orbital (particle or hole) is defined by 
the position of the orbital with respect to the Fermi level 
of the core configuration. It is clear from Fig. 2 that 
the neutron states :/[523]7/2±\ :.[411]l/2±\ z/[413]5/2±^ 
j.[541]l/2±\ ;y[532]5/2±* and u6^' have hole character, 
while i^6t\ 1^63", i/[530]l/2±», i/[402]5/2±», iy[532]3/2±*, 
and z.[514]9/2='=* have particle character. In a simi- 
lar way, the proton orbitals 7t[541]3/2±*, 7r[422]3/2±*, 
7r[301]l/2±*, and 7r[420]l/2±^ and 7r[404]9/2±* can 
be viewed as holes, while 7r[532]5/2±% 7r[411]3/2±', 
7r[541]l/2='=', and 7r[413]5/2='= have particle character (see 
Fig. 3). 



III. RESULTS OF THE ADDITIVITY ANALYSIS 
A. Effective charge quadrupole moments qlo.a 

Table I contains the values of CHF and CRMF effec- 
tive s.p. charge quadrupole moments 920.0 ^ number 
of s.p. orbitals in the vicinity of the deformed shell gaps 
at Z=58 and iV=73 (see Figs. 2 and 3). There is an 
overall excellent agreement between ^lo « values for the 
two mean-field approaches employed. In the majority of 
cases, the uncertainties are small enough to allow deter- 
mination of effective moments to two significant digits. 

The two lowest neutron intruder orbitals 6]"* and 62 * 
show significant signature splitting, and their effective 
charge quadrupole moment values differ by more than 
5%. The extracted values confirm the general expecta- 
tions for the polarization effects exerted by the intruder 
and extruder states [37, 38]. The lowest neutron A/'=6 
orbitals, and 6^*, have 9200 — 0.37 eb, which indi- 
cates that their occupation drives the nucleus towards 
larger prolate deformation. The third intruder orbital, 
6J*, although calculated with relatively poor statistics, 
confirms this trend. 

The proton 7r[404]9/2='=* extruder high-f2 orbitals are 
oblate-driving; they have large negative values of (i'2o,Q- 
Emptying them polarizes the nucleus towards more 
prolate-deformed shapes. Interestingly, their gloa val- 
ues of around — 0.31 eb are close in magnitude to those of 
the M=6 neutron intruders, in line with the findings of 
Ref. [18] that the holes in the proton 59/2 orbitals are as 
important as the particles in the neutron ii3/2 orbitals in 
stabilizing the shape at large deformation. Due to their 
high-il content, the signature splitting of 7r[404]9/2^* 
routhians is extremely small and their ^20 a values are 
practically indistinguishable within error bars. 

Our study indicates that proton ftii/2 states, such as 
7r[541]3/2±* and 7r[532]5/2±' active below and above the 
Z=58 shell gap, respectively, play a significant role in 
the existence of this island of high deformation. Indeed, 
Table I attributes them to effective charge quadrupole 
moments in excess of 0.45 eb - very significant values com- 
pared with other states listed. 

The downsloping orbital 7r[541]l/2~*, originating from 
mixed 7r/7/2 © ^9/2 subshells, carries a large effective 
charge quadrupole moment of more than 0.5 eb. Al- 
though one could expect it to play a role in the formation 
of large prolate deformation, this state appears too high 
in energy (above the Z=58 shell gap) and would there- 
fore always stay unoccupied in most of the configurations 
of interest [19-21]. On the contrary, the strongly prolate- 
driving 7r[550]l/2~* orbital carrying (72o,Q~0-47eb, is al- 
ways occupied in the bands of interest. 

Table I compares the values of (720. « obtained in the 
present study with those from the additivity analysis 
of the SD bands in the ^~150 region [5] based on the 
Skyrme SkP and SkM* energy density functional. Note 
that some of the states, which are of particle character in 
the A~130 region, appear as hole states in the heavier re- 
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TABLE I: Effective s.p. charge quadrupole moments glo.a (in eb) for the s.p. orbitals active in the yl~130 mass region of high- 
and supcrdoformation. Calculations were carried out with CHF+SLy4 and CRMF+NLl approaches. The bare quadrupole 
moments q20,a are also shown for CHF+SLy4. Theoretical errors resulting from the multivariate least-square fit are indicated. 
The results of previous calculations [5] pertaining to the j4~150 mass region are displayed for comparison. 



State 



CHF+SkP CHF+SkM=' 



920, c 



[402] I 

[402] I" 

[411] r 
[411] r 
[411] r 

[411] I " 

[413] r 

[413] I - 



p/h 


y20,a 


p 


-0.44 


p 


-0.44 


h 




h 




h 




h 




h 


-0.16 


h 


-0.13 





CHF+SLy4 




CRMF+NLl 


n „barc 




off 
120,0, 






off 
920, Q 




p 


0.0 


-0 


35±0 


01 


-0 


26±0 


01 


p 


0.0 





34±0 


02 





26±0 


02 


h 


0.0 


-0 


15±0 


02 


-0 


11±0 


02 


h 


0.0 


-0 


12±0 


01 


-0 


06±0 


02 


h 


0.0 


-0 


15±0 


04 





13±0 


03 


h 


0.0 


-0 


11±0 


05 


-0 


12±0 


03 


h 


0.0 


-0 


13±0 


02 


-0 


13±0 


03 


h 


0.0 


-0 


12±0 


03 


-0 


11±0 


02 


h 


0.0 





03±0 


01 





05±0 


01 


h 


0.0 





04±0 


01 





01±0 


02 


P 


0.0 





22±0 


01 





17±0 


01 


P 


0.0 





17±0 


Ul 





19±0 


01 


P 


0.0 





21±0 


03 








P 


u.u 





17±0 


03 









u 
n 


0.0 





19±0 


03 





17±0 


03 


u 
n 


0.0 





24±0 


03 





38±0 


03 


n 


U.U 





35±0 


03 





35±0 


02 


h 


0.0 





37±0 


03 





33±0 


03 


I, 
n 


0.0 





38±0 


01 





40±0 


01 


P 


U.U 





36±0 


01 





36±0 


01 


P 


U.U 





35±0 


05 









h 


-0.08 


u 




05 








P 


-0.13 


-0 


32±0 


01 


-0 


37±0 


01 


P 


-0.13 


-0 


32±0 


01 


-0 


37±0 


01 


P 


0.06 


-0 


05±0 


02 








P 


0.06 





00±0 


01 








P 


0.06 





28±0 


05 








h 


0.20 





33±0 


02 





33±0 


03 


h 


0.22 





34±0 


02 





28±0 


02 


P 


0.28 





43±0 


01 





41±0 


02 


P 


0.36 





56±0 


03 





54±0 


03 


P 


0.40 





58±0 


02 








h 


0.34 





50±0 


01 





48±0 


01 


h 


0.39 





57±0 


01 





50±0 


01 


h 


0.30 





49±0 


05 





47±0 


04 



// [523] 

ly [523] 

v [530] 

u [530] i 

1/ [532] I 

1/ [532] I 
[532] I 

V [532] I 

p [541] i 

u [541] i 



-0.38 
0.38 
-0.18 
-0.15 



h 


0.43 


0.30 


h 


-0.15 


-0.13 


p 


-0.30 


-0.28 


p 


-0.30 


-0.28 


p 


0.11 


0.10 


p 


0.11 


0.10 



TT [301] 

TT [404] I 



[404] I 
[411] I 



TT [411] I 
TT [413] I 
TT [422] I 
TT [422] I 

TT [532] I" 
TT [532] I" 
TT [541] I" 
TT [541] I " 
7r[541]|- 
TT [550] i- 



gion. For these states, conforming to our definitions of co- 
efficients Cq, (Sec. II B), we inverted signs of values shown 
in Table 1 of Ref. [5]. With few exceptions, q values 
are similar in both studies: only for the 7r[301]l/2"'"* or- 
bital does the; difference between A ~ 130 and A ~ 150 
results exceed 0.1 eb. This result strongly suggests that 
the polarization effects caused by occupying/emptying 
specific orbitals are mainly due to the general geometric 



properties of s.p. orbitals and weakly depend on the ac- 
tual parametrization of the Skyrme energy density func- 
tional; minor differences are likely related to interactions 
between close-lying s.p. states. These observations give 
strong reasons for combining the two regions into one, 
and interpreting the entire area of highly and SD rota- 
tional states in the mass range A 128 — 160 within the 
united theoretical framework. 
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The results for <720 a obtained in CHF+SLy4 and 
CRMF+NLl models are indeed very similar (see Table 
I). Only for the zy[402]5/2±* and t/[532]5/2-* orbitals, do 
the differences between glo a values come close to 0.1 eb. 

Table I compares the bare and effective s.p. charge 
quadrupole moments obtained in CHF+SLy4. In the ma- 
jority of cases, these quantities differ drastically, under- 
lying the importance of shape polarization effects. Large 
differences between bare and effective s.p. quadrupole 
moments have also been found in the CHF+SkP and 
CHF-|-SkM* calculations in the A-^ISO region of superde- 
formation [5]. 

B. Effective quadrupole moments and q22,a 

Table II displays the calculated effective s.p. transi- 
tion quadrupole moments gf^, cf. definitions (15-16). 
Based on the additivity principle, these values can be 
used to predict the total charge transition moments Qt {k) 
in highly deformed and SD bands of A ~ 130 nuclei: 

Qt{k) = Qr' + Y,Ca{k)q!%, (30) 

a 

where the calculated CHF+SLy4 value for the core con- 
figuration in ^^^Cs is 

Q™^<' = 7.64eb. (31) 

Since the total calculated values are less precise than 
the relative ones which define the effective s.p. transi- 
tion quadrupole moments one may alternatively use 
in Eq. (30) the measured value [39], 

Q™"P = 7.4(3) eb. (32) 

Theoretical estimates of the total charge transition mo- 
ments Qt{k) allow for predictions of B{E2) values 

B{E2){I ^ I -2,k) = -^e2(/0 20|/-2 0)g2(/c), (33) 

and lifetimes [40]. 

In CHF-|-SLy4, the uncertainties of q'22 a appear to be 
larger than those for glo.a- CRMF-l-NLl, however, 
those uncertainties arc similar. This can be traced back 
to the different 7-softness of potential energy surfaces in 
CHF-FSLy4 and CRMF-l-NLl (see Ref. [18] and refer- 
ences quoted therein for the results obtained in different 
approaches); current analysis revealing large uncertain- 
ties for glf „ suggests that the potential energy surfaces 
are softer (and, thus less localized) in the CHF-|-SLy4 
approach. 

Although the values of g|f „ are generally much smaller 
than a 1 large uncertainties in the determination 
of certain moments 9220 (especially, for i'[411]3/2='=*, 
!/[532]5/2±', vQ^\ 7r[301]'l/2+% and 7r[550]l/2-^ orbitals, 
for which the errors exceed 0.1 cb in the CHF+SLy4 ap- 
proach) can lead to the deterioration of predicted qf^^. 



On the other hand, in many cases the uncertainties in 
(722 a arc smaller than the experimental error bars; hence, 
they are less relevant when comparison with experiment 
is carried out. Currently available experimental data on 
relative transition quadrupole moments agree reasonably 
wen with the CHF-hSLy4 results [20, 21]. 

Table II compares values obtained in CHF-|-SLy4 
and CRMF+NLl models. The results for proton orbitals 
are similar in both approaches: the differences between 
respective qf^^ values do not exceed 0.1 eb. Larger dif- 
ferences are seen for the neutrons: for about 50 per- 
cent of calculated orbitals (i^[402]5/2±^ z/[411]l/2±\ 
^/[411]3/2+^ ^/[413]5/2-^ ;y[530]l/2+*, and i/[532]5/2+*), 
the difference between values in CIIF+SLy4 and 
CRMF+NLl exceeds 0.1 eb. Interaction (mixing) be- 
tween those close-lying states (see Fig. 2), predicted dif- 
ferently in the two approaches, is the most likely reason 
for the deviations seen. 

The results of CHF+SLy4 were compared with exper- 
imental transition moments in Rcfs. [20, 21]. Here, we 
show in Fig. 4 a comparison between CRMF+NLl and 
experiment for the relative transition quadrupole mo- 
ments 5Qi{k) in different highly deformed and SD bands 
in nuclei with Z=57-62 involving «i3/2 neutrons and/or 
gg/2 proton holes. The agreement between experiment 
and theory is quite remarkable with all the experimen- 
tal trends discussed in Refs. [20, 21] well reproduced by 
calculations. One should note that the CRMF and CHF 
results are close to each other. The general pattern of de- 
creasing Qt with increasing Z and N is consistent with 
the general expectation that as one adds particles above 
a deformed shell gap, the deformation-stabilizing effect of 
the gap is diminished. This trend continues until a new 
"magic" deformed number is reached. Such a situation 
occurs when going from ^^^Ce towards Z=62 and A^=80 
(^^^Sm), where a large jump in transition quadrupole mo- 
ment takes place marking the point at which it becomes 
energetically favorable to fill the high-j 7rzi3/2 and vji5/2 
orbitals responsible for the existence of the A 142 SD 
island. 

It is gratifying to see that CRMF+NLl reproduces the 
value of Qt in ^^^Sm based on the ^"^^Ce core (see inset 
in Fig. 4). Earlier on, it was demonstrated in Refs. [9, 
20, 21] that this Qt value can be also reproduced within 
CHF using either a ^^^Ce or a ^^^Dy core. 



C. Effective angular momenta ja (s.p. alignments) 

In this section, we evaluate and interpret the effective 
s.p. contributions to the total angular momentum. Ta- 
ble HI displays effective s.p. angular momenta for 
the s.p. orbitals of interest. The relative uncertainties 
in calculated values arc on average larger than those 
for effective s.p. quadrupole moments. This is due to 
the fact that, on the mean field level, polarization effects 
pertaining to the angular momentum are more complex 
than those for quadrupole moments: they involve not 
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FIG. 4: Experimental (closed symbols with error bars) and calculated (CRMF-l-NLl, open symbols) differential transition 
quadrupole moments for highly deformed bands in Ce, Pr, Nd, Pm, and Sm isotopes. The experimental data were taken from 
Refs. [20, 21] and references quoted therein. The values of SQt for SD band in ^'''^Sm are shown in the inset. Dashed lines are 
drawn to guide the eye. 



only shape changes but also the variations of time-odd 
mean fields [41-43]. For eight proton states calculated in 
both approaches, the mean uncertainties are 0.19fi. and 
0.18ft in CHF+SLy4 and CRMF+NLl, respectively. The 
same holds also for the set of 18 neutron states, where the 
average uncertainties are 0.25ft and 0.20ft in CHF-|-SLy4 
and CRMF+NLl, respectively 

Table III also compares CHF+SLy4 expectation values 

of the s.p. angular momentum j^^"^ = with their 

effective counterparts j^. It is seen that these two quan- 
tities differ considerably. As discussed in Ref . [42] , this is 
due to both shape polarization and time-odd mean-field 
effects. It is also important to remember that, unlike the 
cranked Nilsson scheme, in self-consistent models the ex- 
pectation value of the projection of the s.p. angular mo- 
mentum on the rotation axis j^^'^'^ cannot be extracted 
from the slope of its s.p. routhian versus rotation fre- 
quency [44]. 

Our results indicate that the additivity principle for 
angular momentum alignment does not work as precisely 
as it does for quadrupole moments. This conclusion is 
in line with a similar analysis in the A ~ 60 region of 
superdeformation [45, 46]. A configuration assignment 
based on relative alignments depends on how accurately 
these alignments can be predicted. For example, the 
application of effective (relative) alignment method in 
the A ~ 140 — 150 region of superdeformation requires 
an accuracy in the prediction of relative angular mo- 
menta on the level of ~ 0.3ft and 0.5ft for non-intruder 
and intruder orbitals, respectively [8, 13, 47]. In the 



highly deformed and SD nuclei from the ^ ~ 60 — 80 
mass region, these requirements for accuracy are some- 
what relaxed (see Refs. [48, 49]). We expect that in 
the A ~ 130 region, the relative alignments should be 
predicted with a precision similar to that in the A ^ 
140 — 150 region. However, for a number of orbitals (for 
example, i/[411]3/2±^ z/[532]5/2±', j/6^\ 7r[301]l/2+', 
7r[413]5/2^*, 7r[550]l/2~*), the calculated uncertainties 
in are close to 0.5ft, and this probably prevents reli- 
able assignments based on the additivity principle for the 
configurations involving these orbitals. The situation be- 
comes even more uncertain if several orbitals with high 
uncertainties in are occupied. 

Let us also remark that while theory provides effec- 
tive alignments at a fixed rotational frequency, relative 
alignments extracted from experimental data may show 
appreciable frequency dependence (see for instance Ref. 
[45]). Therefore, for reliable configuration assignments, 
measured relative alignments should be compared with 
calculated ones over a wide frequency range. 



D. Variance and distribution of residuals 

One of the main outcomes of this study is the set of ef- 
fective s.p. moments 520,01 122,0^ 9t,Q; ^"^d ahgnments j^. 
The quality of the additivity principle can be assessed 
by studying the distribution of first moments of resid- 
uals (25), i.e., differences between the self-consistently 
calculated values of physical observables and those ob- 



11 



TABLE II: Effective s.p. charge quadrupole moments glo.a and 922,0 as well as the transition quadrupole moments qt^a (all in 
eb) calculated in CHF+SLy4 and CRMF+NLl. 



State 



CSHF + SLy4 





920, a 






922, a 












920, a 






922, a 






<lt,a 




-0 


35±0 


01 





14±0 


06 


-0 


04±0 


04 


-0 


26±0 


02 


-0 


02±0 


01 


-0 


25±0 


02 





34±0 


02 





08±0 


08 


-0 


38±0 


05 


-0 


26±0 


02 


-0 


07±0 


02 


-0 


22±0 


03 


-0 


15±0 


02 


-0 


24±0 


10 


-0 


01±0 


06 


-0 


11±0 


02 





09±0 


02 


-0 


16±0 


02 


-0 


12±0 


01 





06±0 


06 -0 


16±0 


04 


-0 


06±0 


02 


-0 


17±0 


02 





04±0 


02 


-0 


15±0 


04 





20±0 


20 


-0 


26±0 


12 


-0 


13±0 


03 


-0 


02±0 


03 


-0 


11±0 


03 


-0 


11±0 


05 


-0 


05±0 


24 


-0 


08±0 


15 


-0 


12±0 


03 





02±0 


03 


-0 


12±0 


03 


-0 


13±0 


02 


-0 


05±0 


10 -0 


10±0 


06 


-0 


13±0 


03 


-0 


04±0 


03 


-0 


10±0 


03 





12±0 


03 





12±0 


13 





05±0 


08 





llzbO 


02 





15±0 


03 


-0 


20±0 


03 





03±0 


01 


-0 


00±0 


05 





03±0 


03 





05±0 


01 





00±0 


01 





04±0 


01 





04±0 


01 


-0 


01±0 


05 





05±0 


03 





01±0 


02 


-0 


00±0 


02 





01±0 


02 





22±0 


01 


-0 


21±0 


05 





34±0 


03 





17±0 


01 


-0 


09±0 


01 





22±0 


01 





17±0 


01 


-0 


01±0 


Uo 





18±0 


03 





19±0 


01 





10±0 


01 





13±0 


01 


u 


91 -t-O 


Uo 


n 
u 


91 -4-0 


13 





09±0 


08 























17±0 


03 





03±0 


13 





15±0 


08 























19±0 


03 


-0 


08±0 


20 





24±0 


12 





17±0 


03 


-0 


02±0 


03 





18±0 


03 





24±0 


03 


-0 


01±0 


20 





25±0 


12 





38±0 


03 





00±0 


03 





38±0 


03 





35±0 


03 


-0 


04±0 


13 





38±0 


08 





35±0 


02 


-0 


00±0 


02 





35±0 


03 





37±0 


03 





01±0 


14 





36±0 


08 





33±0 


03 





04±0 


03 





30±0 


03 


Q 


38±0 


01 





21±0 


03 





26±0 


02 





40±0 


01 


Q 


12±0 


01 





33±0 


01 





36±0 


01 





oi±n 


04 





37±0.03 





36±0 


01 





01±0 


01 





37±0 


01 


I) 


:!■")=() 


()."•) 





()()±() 


22 


().:icS±() 

























51±0 


05 


-0 


10±0 


24 





57±0 


14 




















-0 


32±0 


01 





10±0 


04 


-0 


38±0 


02 


-0 


37±0 


01 





02±0 


01 


-0 


38±0 


01 


-0 


32±0 


01 





09±0 


04 


-0 


37±0 


02 


-0 


37±0 


01 





02±0 


01 


-0 


38±0 


01 


-0 


05±0 


02 





10±0 


07 


-0 


10±0 


05 























00±0 


01 





22±0 


07 





12±0 


04 























33±0 


02 


-0 


27±0 


10 





48±0 


06 





33±0 


03 


-0 


13±0 


02 





40±0 


03 





34±0 


02 





14±0 


10 





25±0 


06 





28±0 


02 





16±0 


02 





19±0 


02 





43±0 


01 


-0 


05±0 


05 


-0 


46±0 


03 





41±0 


02 


-0 


04±0 


01 





43±0 


02 





56±0 


03 


-0 


07±0 


09 





60±0 


05 





54±0 


03 





05±0 


03 





51±0 


04 





58±0 


02 


-0 


01±0 


10 





59±0 


06 























50±0 


01 


-0 


05±0 


06 





52±0 


04 





48±0 


01 


-0 


10±0 


01 





54±0 


01 





57±0 


01 


-0 


12±0 


04 





63±0 


03 





50±0 


01 


-0 


10±0 


01 





56±0 


01 





49±0 


05 


-0 


06±0 


22 





52±0 


14 





47±0 


04 


-0 


02±0 


04 





48±0 


04 



CRMF + NLl 



z.[402]| + 



p[402] 
K411]i + 
i/[411]i ' 
^[411]!+' 

^[4ii]r' 

z/[413]|+' 

z.[413l|-' 

u[523]l + 
!/[523]|" 
i^[530]i + 
;.[530]i- 
i'[532]| + 
i^[532]|" 
z^[532]| + 
i/[532]|" 
z.[541]l+ 
i.[541]i 



u 6~ 

per 

V 6., 



7r[301]i+' 

7r[404]| + 
7r[404]|" 
7r[411]| + 
^[411] I 
7r[422]| + 
7r[422]|" 

7r[532]|+ 

7r[532]|" 

7r[541]i" 

7r[541]|+ 

7r[541]| 

7r[550]i-^ 



tained from the additivity principle. For instance, for 
the quadrupole moment Q20-, the quantity of interest is 

AQ20 = - 5Q2o{k). (34) 

a 

Deviations A(522: AQt: and A J arc given by similar ex- 
pressions. Figures 5 and 6 show distributions of these 
deviations. The quality of the additivity principle for 
Q2V1 is shown in the top two panels of Fig. 5. In the CHF 
model, the majority of A(520 values (more than 97.8% of 



the total number) fall comfortably within the interval of 
±0.1 eb. This corresponds to a relative distribution width 
of about 1.3%. In CRMF, the distribution is even nar- 
rower, with more than 90% of AQ20 values falling within 
the ±0.05 eb interval, or less than 0.7% of the total value. 

The results for the total angular momentum are shown 
in the bottom panels of Fig. 5. In CRMF, the distribu- 
tion of deviations is very narrow, with only 10% of the 
cases differing by more than ±ft/2. The CHF histogram 
is somewhat wider, but more than 90% of deviations fall 
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TABLE III: Effective s.p. angular momentum alignments j'^ 
(in h) of the active orbitals calculated in CHF+SLy4 and 
CRMF+NLl. In the second column, the bare s.p. angular 
momenta j^^'^^, calculated with CHF+SLy4 are also shown. 



State 
[Afn.A]Q' 



U 

I 2 

ir 
1+' 

2 



V [402] I +' 

V [402] I ~ 

^[411]^' 

V [411] i 

u [411] §+■ 

r ..... T — 

V 

V [413] 

V [413] I 

[523] 

V [523] I" 
[530]!"^ 

V [530] i"' 
[532] I 
[532] I" 

V [532] I +■ 
[532] I" 

V [541] i+' 
[541] i- 



[301] 



1 + 



TT [404] I + 
TT [404] I" 

^ [411] r 

^ [411] I 



TT 



[413] I 
[422] I + 



TT 

TT [422] 
TT 



TT 



[532] I + 
[532]|- 

[541] r 

[541] 1+ 
[541] I - 
[550] i 



CHF+SLy4 

■bare 
Jet J a 



-0.528 
-0.493 
0.411 
0.380 
-0.092 
0.077 
-0.316 
-0.428 

-0.908 
-0.974 
1.548 
0.564 
0.171 
0.835 
-0.331 
0.417 
1.793 
0.466 

4.840 
4.031 
2.662 



0.58±0.14 
0.51±0.20 
0.67±0.24 
0.40±0.15 
1.72±0.46 
0.56±0.57 
-0.10±0.23 
0.12±0.30 

-l.lOiO.lO 
-1.19±0.12 
1.19±0.11 
0.88±0.11 
-0.34±0.30 
0.44±0.31 
-0.89±0.46 
-1.06±0.46 
0.92±0.31 
0.89±0.32 

4.78±0.08 
3.42±0.11 
0.77±0.50 



-0.432 
-0.719 
-0.719 
-0.249 
-0.062 
-0.539 
-0.315 

0.510 
-0.253 
-0.022 

0.944 

1.743 
-0.057 

2.819 



1.23±0.55 
-0.00±0.09 
-0.00±0.09 
0.81±0.18 
0.65±0.16 
-1.52±0.53 
-0.19±0.25 
-0.84±0.23 
-0.90±0.13 
-0.67±0.20 
1.75±0.23 
1.57±0.13 
-0.54±0.10 
2.99±0.52 



CRMF-hNLl 

•eff 
Jo. 



0.47±0.15 
0.38±0.26 
0.64±0.17 
0.09±0.16 
1.35±0.29 
1.08±0.29 
0.44±0.27 
0.14±0.26 

-1.24±0.12 



-0.95±0.29 
-1.29±0.30 
0.95±0.25 
-0.34±0.28 

4.59±0.08 
3.15±0.10 



0.09±0.09 
0.11±0.09 



-0.21±0.27 
-0.38±0.24 
-1.11±0.16 



1.18±0.11 
-0.48±0.11 
2.86±0.40 



20 



10 



20 
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183 bands 



n r* \ n n n I n n n 
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105 bands 
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FIG. 5: Histogram of differences between self-consistent val- 
ues obtained in CHF and CRMF and those given by the ad- 
ditivity formula [see e.g. Eq. (34)]. The results for Q20 are 
shown in the two upper panels and those for the total angular 
momentum are displayed in the two lower panels. 
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FIG. 6: Similar to Fig. 5 except for Q22. 



within the ±fi,/2 interval. Taking into consideration that 
the experimental spins of highly deformed and SD bands 
are often assigned with uncertainties that are multiples 
of h, our results give considerable encouragement for the- 
oretical interpretations based on the method of relative 
(effective) alignments [8, 13, 47]. 

In CHF and CRMF, the distributions of deviations of 
charge quadrupole moments Q22 (Fig. 6) are relatively 
narrow. Again, for CRMF, nearly 95% of deviations fall 
within ±0.025 eb, and 98% fall within ±0.1 eb. For CHF, 



the distribution of deviations is somewhat wider, with 
more than 90% of deviations falling within the ±0.2 eb 
interval. 

We interpret these results as a strong indication that 
the additivity principle works fairly well in self-consistent 
cranked theories. While distributions of deviations in 
Q20 and Q22 are rather similar in the CHF+SLy4 and 
CRMF+NLl models (see top of Fig. 5 and Fig. 6), de- 
viations in angular momentum differ between these two 
approaches. Considering that (i) the uncertainties in j"^ 
are similar in both methods (Sec. IIIC), and (ii) shape 
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polarization effects are not that different (Sec. IIIB), one 
can conclude that the observed difference is due to the 
polarization of time-odd mean fields. However, the de- 
tailed investigation of this effect is beyond the scope of 
this study. 

IV. CONCLUSIONS 

The additivity principle in highly and SD rotational 
bands of the ^ ~ 130 mass region has been studied within 
the cranked Hartree+Fock theory based on the SLy4 en- 
ergy density functional and the cranked relativistic mean- 
field theory with the NLl Lagrangian. The main results 
can be summarized as follows: 

• The two sets of effective s.p. charge quadrupolc mo- 
ments ^22 > transition quadrupole moments 
gf*, and effective angular momenta j'^^ have been 
produced. This rich output allows for an easy and 
simple determination of transition quadrupole mo- 
ments Qt in highly deformed and SD bands in the 
A ^ 130 mass region. In some cases, configuration 
assignments based on the relative (effective) align- 
ment method can be done based on the calculated 
values of effective angular momenta j^^ (see, how- 
ever, Sec. IIIC). 

• Our statistical analysis of distributions of residuals 

confirms that the additivity principle is well ful- 
filled in the self-consistent approaches that properly 
take into account polarization effects. 

• The contribution from the triaxial degree of free- 
dom to the transition rjuadrupole moment is usu- 
ally small, but it cannot be ignored when aiming at 
a quantitative reproduction of experimental data. 
The average magnitude of values is greater in 
the CHF-hSLy4 model than in CRHF-^NLl, thus 
suggesting that the potential energy surfaces pro- 
duced in the former model are more 7-soft. 



• For the majority of s.p. orbitals, there is a consid- 
erable difference between the effective and bare ex- 
pectation values of one-body operators. This indi- 
cates the importance of polarization effects (shape 
polarization for quadrupole moments and the shape 
and time-odd-mean-field polarization for angular 
momentum alignment). 

• With very few exceptions, there is a great deal of 
consistency between CHF and CRMF results for 
the effective s.p. moments and alignments. 

So far, the additivity principle has been investigated only 
for highly deformed or SD bands in the yl^ 130- 150 mass 
region. It would be interesting to extend such studies to 
other high spin structures. The most promising candi- 
dates arc: (i) terminating bands in the ^ ^ 110 mass 
region characterized by very weak pairing and apprecia- 
ble 7-softness [14], and (ii) SD rotational bands in the 
^ ^ 60 and ^ ^ 80 mass regions of superdeformation 
[46]. Work along these lines is in progress. 



V. ACKNOWLEDGEMENTS 

Our study was inspired by the experimental work of 
Laird et al. [21]. Stimulating discussions with Mark 
Riley are gratefully acknowledged. The work was sup- 
ported in part by the U.S. Department of Energy 
under Contract Nos. DE-FG02-96ER40963 (University 
of Tennessee), DE-AC05-00OR22725 with UT-Battelle, 
LLC (Oak Ridge National Laboratory), and DE-FG05- 
87ER40361 (Joint Institute for Heavy Ion Research), and 
DE-FG02-07ER41459 (Mississippi State University); by 
the Latvian Scientific Council (grant No. 05.1724); by 
the Polish Ministry of Science; by the Academy of Fin- 
land and University of Jyvaskyla within the FIDIPRO 
programme; and by the European Union Social Fund and 
the research program Pythagoras II - EPEAEK II, under 
project 80861. 



[1] I. Ragnarsson, Phys. Lett. B 264, 5 (1991). 

[2] G. de Prance, C. Baktash, B. Haas, and W. Nazarewicz, 

Phys. Rev. C 53, R1070 (1996). 
[3] L. B. Karlsson, I. Ragnarsson, and S. Aberg, Phys. Lett. 

B 416, 16 (1998). 
[4] J. Zhang, Y. Sun, L. L. Riedinger, and M. Guidry, Phys. 

Rev. C 58, 868 (1998). 
[5] W. Satula, J. Dobaczowski, J. Dudck, and W. 

Nazarewicz, Phys. Rev. Lett. 77 5182 (1996). 
[6] S. Aberg, L.-O. Jonsson, L. B. Karlsson, and I. Ragnars- 
son, Z. Phys. A 358, 268 (1997). 
[7] L. B. Karlsson, I. Ragnarsson, and S. Aberg, Nucl. Phys. 

A 639, 654 (1998). 
[8] A. V. Afanasjev, G. Lalazissis, and P. Ring, Nucl. Phys. 

A 634, 395 (1998) 



[9] G. Hackman, R. V. F. Janssens, E. F. Moore, D. Nisius, 
1. Ahmad, M. P. Carpenter, S. M. Fischer, T. L. Khoo, T. 
Lauritson, and P. Rcitcr, Phys. Lett. B 416, 268 (1998). 

[10] S. T. Clark, G. Hackman, R. V. F. Janssens, R. M. Clark, 
P. Fallon, S. N. Floor, G. J. Lane, A. O. Macchiavclh, J. 
Norris, S. J. Sanders, and C. E. Svensson, Phys. Rev. 
Lett. 87, 172503 (2001). 

[11] W. Nazarewicz, R. Wyss, and A. Johnson, Nucl. Phys. 
A 503, 285 (1989). 

[12] A. V. Afanasjev, J. Konig, and P. Ring, Nucl. Phys. A 
608, 107 (1996). 

[13] I. Ragnarsson, Nucl. Phys. A 557, 167c (1993). 

[14] A. V. Afanasjev, D. B. Fossan, G. J. Lane, and 1. Rag- 
narsson, Phys. Rep. 322, 1 (1999). 

[15] W. Satula and R. Wyss, Rep. Prog. Phys. 68, 131 (2005). 



14 



[16] G. Stoitcheva, W. Satula, W. Nazarowicz, D.J. Dean, M. 
Zalcwski, and H. Zdunczuk, Phys. Rev. C 73, 061304(R) 
(2006). 

[17] R. Wyss, J. Nyberg, A. Johnson, R. Bengtsson, and W. 

Nazarewicz, Phys. Lett. B 215, 211 (1988). 
[18] A. V. Afanasjev and I. Ragnarsson, Nucl. Phys. A 608, 

176 (1996). 

[19] E. S. Paul, P. T. W. Choy, C. Andreoiu, A. J. Boston, A. 
O. Evans, C. Fox, S. Gros, P. J. Nolan, G. Rainovski, J. 
A. Sampson, H. C. Scraggs, A. Walker, D. E. Appelbe, 

D. T. Joss, J. Simpson, J. Gizon, A. Astier, N. Buforn, 

A. Prevost, N. Redon, O. Stezowski, B. M. Nyako, D. 
Sohler, J. Timar, L. Zolnai, D. Bazzacco, S. Lunardi, C. 
M. Petrache, P. Bednarszyk, D. Curien, N. Kintz, and I. 
Ragnarsson, Phys. Rev. C 71, 054309 (2005). 

[20] M. A. Riley, R. W. Laird, F. G. Kondev, D. J. Hartley, D. 

E. Archer, T. B. Brown, R. M. Clark, M. Devlin, P. Fal- 
lon, I. M. Hibbert, D. T. Joss, D. R. LaFosse, P. J. Nolan, 
N. J. O'Brien, E. S. Paul, J. Pfohl, D. G. Sarantites, R. 
K. Sheline, S. L. Shepherd, J. Simpson, R. Wadsworth, 
M. T. Matov, A. V. Afanasjev, J. Dobaczewski, G. A. 
Lalazissis, W. Nazarewicz, and W. Satula, Acta Phys. 
Polonica B 32, 2683 (2001). 

[21] R. W. Laird, F. G. Kondev, M. A. Riley, D. E. Archer, T. 

B. Brown, R. M. Clark, M. Devlin, P. Fallon, D. J. Hart- 
ley, I. M. Hibbert, D. T. Joss, D. R. LaFosse, P. J. Nolan, 
N. J. O'Brien, E. S. Paul, J. Pfohl, D. G. Sarantites, R. 
K. Sheline, S. L. Shepherd, J. Simpson, R. Wadsworth, 
M. T. Matev, A. V. Afanasjev, J. Dobaczewski, G. A. 
Lalazissis, W. Nazarewicz, and W. Satula, Phys. Rev. 
Lett. 88, 152501 (2002). 

[22] J. Dobaczewski and P. Olbratowski, Comp. Phys. Comm. 

158, 158 (2004). 
[23] P. Ring, A. Hayashi, K. Kara, H. Emling, and E. Grosse, 

Phys. Lett. HOB, 423 (1982). 
[24] I. Hamamoto and Zheng Xing, Phys. Scr. 33, 210 (1986). 
[25] D.A Varshalovitch, A.N. Moskalev, and V.K. Kersonskii, 

Quantum Theory of Angular Momentum (World Scien- 
tific, Singapore, 1988). 
[26] J. Dobaczewski and J. Dudek, Comput. Phys. Commun. 

102, 166 (1997). 
[27] J. Dobaczewski and J. Dudek, Comput. Phys. Commun. 

102, 183 (1997). 
[28] D. R. Inghs, Phys. Rev. 96, 1059 (1954). 
[29] P. R. Bevingtori, Data Reduction and Error Analysis 

for the Physical Sciences, (McGraw-Hill Publ. Co., New 

York, 1969). 

[30] C. H. Lawson and R. J. Hanson, Solving Least Squares 



Problems, (Prentice-Hall, New York, 1974). 
[31] E. Chabanat, P. Bonche, P. Haensel, J. Meyer, and F. 

Schaeffer, Nucl. Phys. A 635, 231 (1998). 
[32] Mladen Matev, Self-Consistent Description of Rotational 

Properties of Highly Deformed States in Atomic Nuclei 

Far From Stability, PhD thesis. University of Tennessee, 

2003. 

[33] W. Koepf and P. Ring, Nucl. Phys. A 493, 61 (1989) 
[34] J. Konig, and P. Ring, Phys. Rev. Lett. 71, 3079 (1993). 
[35] P.-G. Reinhard, M. Rufa, J. Maruhn, W. Greiner, and J. 

Friedrich, Z. Phys. A 323, 13 (1986) 
[36] A. V. Afanasjev, T. L. Khoo, S. Frauendorf, G. A. 

Lalazissis, and I. Ahmad, Phys. Rev. C 67, 024309 

(2003). 

[37] S. Frauendorf and F. R. May, Phys. Lett. B 125, 245 
(1983). 

[38] Y.S. Chen, S. Frauendorf, and G.A. Leander, Phys. Rev. 
C 28, 2437 (1983). 

[39] R.M. Clark, I.Y. Lee, P. Fallon, D.T. Joss, S.J. Aszta- 
los, J. A. Becker, L. Bernstein, B. Cederwall, M.A. Dele- 
planque, R.M. Diamond, L.P. Farris, K. Hauschild, W.H. 
Kelly, A.O. Macchiavelh, P.J. Nolan, N. O'Brien, A.T. 
Semple, F.S. Stephens, and R. Wadsworth, Phys. Rev. 
Lett. 76, 3510 (1996). 

[40] W. Nazarewicz and I. Ragnarsson, Nuclear Deforma- 
tions, in Handbook of Nuclear Properties, ed. by D.N. 
Poenaru and W. Greiner, (Clarendon Press, Oxford), 
1996, p. 80. 

[41] J. Dobaczewski and J. Dudek, Phys. Rev. C 52, 1827 
(1995). 

[42] A. V. Afanasjev and P. Ring, Phys. Rev. C 62, 031302(R) 
(2000). 

[43] D. Vretenar, A. V. Afanasjev, G. A. Lalazissis, and P. 
Ring, Phys. Rep. 409, 101 (2005). 

[44] B. GaU, P. Bonche, J. Dobaczewski, H. Flocard, and P.- 
H. Heenen, Z. Phys. A 348, 189 (1994). 

[45] J. Dobaczewski, Nuclear Structure'98, AIP Conf. Proc. 
481, edited by C. Baktash (American Institute of 
Physics, New York, 1999), p. 315; rmcl-th/9811043. 

[46] J. Dobaczewski, W. Satula, W. Nazarewicz, and C. Bak- 
tash, to be published. 

[47] C. Baktash, B. Haas, and W. Nazarewicz, Ann. Rev. 
Nucl. Part. Sc. 45, 485 (1995). 

[48] A. V. Afanasjev, I. Ragnarsson, and P. Ring, Phys. Rev. 
C 59, 3166 (1999). 

[49] A. V. Afanasjev and S. Frauendorf, Phys. Rev. C 71, 
064318 (2005). 



